Electronic Structure of gated graphene and graphene ribbons 
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We study the electronic structure of gated graphene sheets. We consider both infinite graphene 
and finite width ribbons. The effect of Coulomb interactions between the electrically injected carriers 
and the coupling to the external gate are computed self-consistently in the Hartree approximation. 
We compute the average density of extra carriers, n2D, the number of occupied subbands and 
the density profiles as a function of the gate potential Vg. We discuss quantum corrections to 
the classical capacitance and we calculate the threshold Vg above which semiconducting armchair 
ribbons conduct. We find that the ideal conductance of perfectly transmitting wide ribbons is 
proportional to the square root of the gate voltage. 
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I. INTRODUCTION 

One important ingredient of semiconductor micro- 
electronics is the electrical tunability of the resistance 
through the field effect. The recent demonstration of field 
effect in graphene^iii^ has opened a new research venue. 
Graphene is different from conventional semiconductors 
for a number of reasons. First, graphene is a truly two 
dimensional atomically thin layer of carbon atoms. Sec- 
ond, neutral graphene is a semimetal with zero density of 
states at the Fermi energy and zero gap, the two scales 
that shape the properties of metals and semiconductors. 
Third, the electronic structure close to the Fermi energy 
has a conical shape with perfect electron-hole symmetry 
and an internal valley symmetry, isomorphic to that of 
two dimensional massless Dirac fermions®. As a result, 
most of the standard lore on electronic and transport 
properties of low dimensional semiconductors needs to 
be revisited^i^i^iiSi. 

The presence of a band-gap in the electronic struc- 
ture is crucial in the design of low dimensional struc- 
tures and in the achievement of large on-off resistance 
ratios in field effect transistors (FETs). However, even 
at very low temperatures, two dimensional graphene 
shows a rather low resistance at the charge neutral- 
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at which the density of states at the 



Fermi energy is vanishingly small. A number of phys- 
ical mechanisms that result in a gap in graphene based 
systems have been proposed: dimensional confinement 
into the so called graphene ribbonsiiii^ii^ii^i^, interlayer 
coupling^®'^ - and spin polarizatior^i^. In this paper we fo- 
cus on graphene ribbons, stripes of graphene with finite 
width W. Ideal graphene ribbons with edges along the 
crystallographic axis fall into two categories: zig-zag and 
armchair—. Only the former can present a gap depend- 
ing on their width W. Recent transport experiments with 
graphene ribbonsi^ confirm that ribbons with W = 2Qnm 
present a thermally activated conductivity that indicates 
the presence of a gap. Most likely, imperfections on the 
edges of real ribbons make it hard to fabricate ribbons 
that are intrinsecally metallic. 



In this paper we explore the relation between gate volt- 
age Vq and injected density in graphene-based field effect 
devices shown in figure 1. A graphene layer, of width W, 
lies on top of an insulating slab of thickness d which lies 
above the metallic gate. Application of a gate voltage 
Vg injects carriers in the graphene layer accompanied by 
a corresponding change in the metal. The main results 
of the paper are the following: (i) we show that Vq and 
the average two dimensional (2D) carrier density in the 
system, n2D , satisfy the relation: 



Vg = 7^en2D 



VQ{n2D) 



(1) 



where the first term arises from the classical electro- 
static interaction and the second term, VQ{n2D), arises 
from quantum mechanical effects associated to the band 
structure of graphene. The classical contribution is dom- 
inant but it depends only on the geometry of the system. 
Although smaller, the second term contains informa- 
tion specific to the electronic structure of the systen>2S. 
(ii) we find that in 2D graphene the quantum contri- 
bution scales like Vq cx y/n2D, and from its measure- 
ment it would be simple to extract the linear slope of the 
graphene bands. (Hi) in the case of ribbons we calculate 
the VG(n2D) curves for different widths W and we find 
the threshold gate above which semiconducting ribbons 
become metallic, (iv) We argue that the ideal conduc- 
tance, G, in wide ribbons scales as G = \/Vg 

The rest of this manuscript is organized as follows. In 
section II we describe the Hartree formalism that we use 
to calculate the electronic structure of gated graphene de- 
vices as well as the link between gate voltage and chem- 
ical potential. In section HI we present the simple an- 
alytical solution of the self-consistent Hartree equation 
for planar graphene {W — oo). We obtain an expression 
for the quantum capacitance of graphene and discuss the 
experimental implications. In section IV we present the 
numerical solution of the self-consistent Hartree equation 
for finite width armchair ribbons, both semiconducting 
and metallic, and we obtain n2D and the conductance G 
as a function of W, and Vg- In section V we discuss the 
implications and present our conclusions. 
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FIG. 1: (Color online), (a) Graphene based FET device: 
lateral view. A graphene ribbon of width W is separated 
from the metallic gate by an insulating slab of thickness d. 
The real charges of the metal are accumulated in the metal- 
insulator interface underneath the ribbon. Image charges lie 
further down the surface, at a distance 2d from the graphene 
(b) Top view of an armchair graphene ribbon. The system is 
infinite along the vertical direction {y axis) . (c) Detail of the 
super-unit cell defining the periodic one dimensional armchair 
ribbon. 



II. THEORETICAL FRAMEWORK 

A. Hamiltonian 

We write the Hamiltonian for electrons in graphene as 
the sum of two parts: Hq + U. The first term describes 
the electronic structure of neutral graphene. We ap- 
proximate Hq by the standard one-orbital tight-binding 
approximation for both the 2D graphene layer— and 
graphene ribbon o^^i^^'^'^ . The model is completely de- 
fined by the positions of the atoms, the first-neighbours 
hopping, t = —2.5eV, the lattice constant a = 2A2A, 
and the average density. The lattice of two dimensional 
graphene is characterized by a crystal structure defined 
by the two unit vectors a± = | (v^, ±l) with a two atom 
basis. We assume that the graphene atoms lie in the 
z = plane. In the case of neutral graphene, the average 
density is one electron per orbital that corresponds to a 
Fermi energy Ep = 0, since we take the on-site energy 
level equal to zero. The armchair ribbons considered in 
this paper are generated by one dimensional repetition of 
the super-cell of fig. lb, along the y axis, with periodicity 
V^a. The supercell itself is constructed as the repetition 
of along the direction (1, 0) of a block with 4 atoms (fig- 
ure Ic). If the supercell has N blocks, the width of the 
ribbon is given hy W = Na and the number of atoms in 
the unit cell as AN. As a rule of thumb, it is useful to 



write W ~ (iV"/4)nm. 

The second term U describes the Coulomb interaction 
between the extra charges in the system: 

U = eJ 5h{r)c^,.,t{T)dr + |^ j ^-!&I!^drd? (2) 

where 5h{r) is the operator describing the departure of 
the local the electronic density from charge neutrality: 

<5n(r)^^|0,(^)P(cLc/,.-"o) • (3) 

Here the sum runs over all the atoms in the lattice and 
0/(r) denotes the tTz atomic orbital centered around the 
atom / and cj ^ is the second quantization operator that 
creates one electron with spin a in the orbital 4>i{r), lo- 
calized around atom /. In Eq[2]no = 1/2 is the number 
of electrons per site and per spin in neutral graphene. 

The electrostatic potential created by the extra charges 
in the metallic gate is denoted by 4>ext{r)- Considered 
as a whole, the graphene and the metal layer form a 
neutral system. Therefore, the extra carriers in one side 
are missing in the other. Since the charges in the metallic 
gate move as to cancel the electric field inside the metal, 
they depend on the density distribution in the graphene 
electrode, which in turn depends on (pexti'f^- This self- 
consistent problem is solved using the image method: for 
a given charge density profile in the graphene, e{5h{r)) 
the potential created by the corresponding extra charges 
in the metal is given by a distribution of fictitious image 
charges inside the metal: —e5nim{r) = —eSn{f—2d) with 
d= (0,0, d). 

V..,ir)^-U-^^^dP . (4) 

With this method the potential evaluated at the metal 
gate (z = 0) is exactly zero. Therefore, the potential 
difference between the metal-insulator interface and the 
graphene layer is the potential evaluated at the graphene 
layer. 



B. Hartree approximation 

The electronic repulsion is treated in the Hartree ap- 
proximation. Therefore, the electrons feel the electro- 
static potentials created by both the gate and themselves: 

U^Usc^eJ 6n{f) (14.t(r) + Vsc(r)) df (5) 

The argument of the integral, Vr = Vext + Vsc, can be 
written as 

Vrir) = - f {Sn{r'))ICi\f-r'\,d)df' (6) 
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where 



2d\ 



(7) 



In the Hartree approximation Hq + Use is a one-body 
Hamiltonian that can be represented in the basis of local- 
ized TT^ atomic orbitals. The _ffo part yields the standard 
single orbital tight-binding Hamiltonian. In contrast, the 
Use introduces site-dependent shifts in the diagonal ma- 
trix elements: 



Ui ^ eVi ^-'^ qjvij 



(8) 



where qj is the average excess charge in site J and 
vjj = J drd'F\(j)i{r)\'^\(j)j{r')\'^IC{\r — T^l^d), is the poten- 
tial created by an unit charge located at site J on site /. 
At this point we adopt an approximation which permits 
to evaluate Vjj without a detailed model of 0/ (r ) . Since 
the atomic orbitals are highly localized, for \Ri — Rj\ > a 
we can approximate vij by the first term in the multi- 



polar expansion: vjj 



We have 



\R,-Rj\ \R,-Rj-2d\ 

verified that the potential in a given site is quite inde- 
pendent of the approximation adopted to evaluate the 
I = J contribution, for which the monopolar approxi- 
mation fails. Without sacrificing accuracy, we adopt the 
simplest strategy of removing that term from the sum. 
The Schrodinger equation reads 



iHo + Usc)\a) 



(9) 



where the eigenvectors \a) are linear combinations of 
atomic orbitals: |a) = ^jC/_q|(^/,(7), and the coeffi- 
cients are independent of the spin because we only con- 
sider non-magnetic solutions. In turn, the distribution of 
excess charges depends on the eigenvalues and eigenvec- 
tors through the equation 



(10) 



where /(eq) is the Fermi function. 

Equations ^ and ^T0\ need to be solved itera- 
tively. In each step of the iterative procedure we have an 
input distribution of extra charge, qi which results in a 
potential (eq. ([8])), which defines a Hamiltonian whose 
eigenstates (eq. ([9])) result in an output distribution of 
extra charge (eq. (fTO]) . The converged solution is such 
that the input and output distribution of charges are the 
same. We refer to that distribution and the correspond- 
ing Hamiltonian as the self-consistent solution. 



C. Gate voltage, chemical potential and quantum 
capacitance 

In real field effect devices, application of a gate volt- 
age Vg results in a change of carrier density in the active 



layer. In our theoretical framework, we define eVc as 
the the chemical potential difference between the metal- 
lic backgate and the graphene ribbon, necessary to ac- 
commodate extra carriers (either holes or electrons) in 
graphene and remove them from the metallic backgate: 



eVG{n2D) = f-graphe 



t^Tnetal 



ill) 



where n2D is the average two-dimensional density of ex- 
tra carriers. In general, there are two contributions to the 
dependence of the chemical potential on n2D ■ One one 
side, the presence of the electrostatic potential created by 
extra carriers present both in graphene and the metallic 
backgate shifts and modifies the bands. We refer to the 
new bands as the Hartree bands. On the other, addition 
of new carriers involves an additional shift of the chemi- 
cal potential with respect to the Hartree bands. Whereas 
the first contribution has a classical origin, the second is 
a consequence of the Pauli principle and we refer to it 
as quantum contribution to the capacitance. We neglect 
a third type of contribution , arising from the modifica- 
tion of the bands due to density dependent exchange and 
correlation contributions^. 

The extra carriers in the metallic gate are concentrated 
on the surface so that the hulk chemical potential and the 
hulk energy bands are shifted by the same amount: there 
is no quantum contribution to the capacitance coming 
from the bulk metal. Things are different in the graphene 
layer. In this case, it makes no sense to talk about bulk 
and surface as different objects. In graphene, the chemi- 
cal potential has to shift relative to the modified bands in 
order to accommodate extra carriers, resulting in a non- 
zero quantum contribution to the capacitance. Without 
loss of generality we take the metallic side as the ref- 
erence for the electrostatic contribution to the chemical 
potential. This permits to write 



eVa — f-graphene{n2D) — f^graphene{0) = Sfi 



(12) 



From an operational point of view, we consider the sys- 
tem to be at zero temperature, and we take the chemi- 
cal potential of graphene at a given density as the lowest 
unoccupied eigenstate of the self-consistent Hamiltonian. 
Identical results are obtained from the equation 



n2D = 



Psc{E)dE 



(13) 



that relates the chemical potential, the extra density via 
the density of states of the self-consistent Hamiltonian, 
Psc{E). 



III. CLASSICAL AND QUANTUM 
CAPACITANCE OF PLANAR GRAPHENE 

In the case of 2D graphene all the atoms of the lattice 
are equivalent and therefore the extra charge qi and the 
effective potential Uj are independent of the location. As 



a matter of fact, for an average extra density n2_D, the 
electrostatic potential is Uj = ^^^n2Dd- Therefore, the 
Hartree bands of charged graphene are equal to the bare 
bands plus a rigid shift: 



13.95 - 



E,[k) = el{k) + n2Dd 



(14) 



In figure 2 (left panels) we plot both the bare bands 
e°(A;) and the Hartree bands E^{k) for d = 300nm and a 
n2D = lO^^cm"^. We see that Hartree bands are shifted 
upwards. Accordingly, the shift in the chemical poten- 
tial needed to accommodate this extra charge will have 
to account for the electrostatic shift of the bands. On 
top of that, the chemical potential needs to move away 
from the Dirac point, which is also shifted electrostati- 
cally. In order to calculate this second contribution, we 
restrict ourselves to the energy region around Ep = 
in which the density of states is linear (linear bands in 
2D). This is a very good approximation since we consider 
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FIG. 2: (Color online), (a): Bare (below) and Hartree 
(above) bands, close to the Dirac point, of two dimensional 



n2D « n° ~ 3.9 IO^^cto ^, the equilibrium electronic graphene with nao = lO^^cm-^ The vertical lines represent 



density in graphene. Therefore we write the density of 
states of the self-consistent Hamiltonian: 



Psc{E) ^ po[E n2Dd 



(15) 



Wxth: po{E) = j^\E\ 

After some simple steps we write the equation relating 
the change in chemical potential dfi 



n2D = / psc{E)dE 
Jo 



which yields 



1 



n2D 



Sp, n2Dd 

e 



(16) 



(17) 



From here we obtain one of the important results of this 
work: 



eVa = dn2D + \t\\ = -r; h eVq (18) 

Thus, we write the gate voltage as the sum of two 
terms. The first is the standard electrostatic contribution 
whereas the second is related to the density of states of 
graphene. We refer to the second, Vq, as the the quan- 
tum contribution to the gate voltage (or inverse capaci- 
tance). In most of devices so far^'^''*'^*' d ~300 nm , which 
makes the electrostatic contribution much larger than the 
quantum contribution. In figure 2b we plot Vg vs n2D, 
according to equation . In spite of the last term of e- 
(HH]), the curve looks like a straight line. As in refs. HF 
we have done a linear fit of fig. 2b obtaining n2D = QfVc 



the chemical potential of the charged system, eVa, which is 
the sum of two contributions, eVc, the rigid shift of the bare 
bands to yield the Hartree bands, and eVg, the shift of the 
chemical potential with respect to the Hartree bands, the 
quantum contribution to the gate. Dotted horizontal lines 
represent the bare and charged Dirac points and dashed hor- 
izontal lines the bare and charged chemical potential, (b) 
Vq = Vc + Vq versus density for 2D graphene in a FET with 
d =300 nm and e = 3.9. (c) Va = Vc + Vq (solid line) and 
Vc (dashed line) versus density for 2D graphene in a FET 
with d =30 nm and e = 47 (see for instance |22) (d) Quantum 
contribution to the gate, Vq versus 7120. 



a thinner dielectric {d =30 nm) with a much larger di- 
electric constant e = 47—. We see how in this case the 
gate voltage (solid line) is signifcantly different from the 
electrostatic contribution (dashed line). In figure 2d we 
plot the quantum contribution, which is independent of 
e and d, alone. Although much smaller than the classi- 
cal term, it can be larger than 0.2 Volt and it should be 
possible to measure it. Therefore, the independent mea- 
surement of n2D, via classical Hall effect, and Vq could 
provide a direct measurement of the slope of the bands 
in the linear region if the experimental results are fitted 
using eq. p8)) . 



with a = 0.718 10"l^-i 



in very good agreement 



with the experiments. In order to reduce the electrostatic 
contribution as much as possible, in figure 2c we consider 



IV. SELF-CONSISTENT ELECTRONIC 
STRUCTURE OF ARM-CHAIR RIBBONS 

A. Bands and charge profiles 

We now consider finite width armchair ribbons. They 
are different from two dimensional graphene both be- 
cause their electronic structure and their electrical capac- 
itance. Armchair ribbons are metallic (figure 3a, c) when 
the width of the sample has the form W=(3M + l)ao. 
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with M an integer, and insulating otherwise (figure 3b, d). 
In both case the electronic structure is different of the 
semimetallic behaviour of 2D graphene. On the other 
side, the presence of edges breaks translational invariance 
along the direction perpendicular to the ribbon axis. The 
edges provide a natural surface to accommodate the ex- 
tra electrons. Therefore, the local potential and the extra 
charge profile are expected to vary as a function of the 
distance to the ribbon edge, keeping their translational 
invariance along the ribbon axis. This makes it harder 
to separate classical and quantum contributions to the 
capacitance in simple terms. 

In figure 3 we show both the self-consistent (upper pan- 
els) and the neutral (lower panels) bands for electrically 
doped ribbons in two cases, semiconducting and a metal- 
lic (iV = 9 and iV = 10 respectively). Notice that charge 
neutral bands have electron-hole symmetry, like the semi- 
conducting carbon nanotubes^^. In those figures we plot 
the self-consistent chemical potential in the upper panels 
and the naive Fermi energy obtained upon integration 
of the neutral density of states. As in the case of 2D 
graphene, the electrical injection of extra carrier results 
in a large shift of the bands, due to the electrostatic inter- 
actions. In contrast to the 2D case, the inhomogeneity of 
the electronic density and the electrostatic potential re- 
sult in a moderate change of the the shape of the bands. 
Notice for instance that for the = 9 ribbon the chem- 
ical potential intersects 2 self-consistent bands whereas 
the naive Fermi energy intersects 3 neutral bands. 

As in the two dimensional case, the self-consistent 
chemical potential can be though as the sum of a large 
electrostatic shift and a smaller quantum shift with re- 
s' 
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FIG. 3: (Color online). Energy bands for A'^ = 9 (left panels) 
and A^ = 10 (right panels). The Hartree (bare) bands are 
shown above (below) . The corresponding densities are n2D ~ 
7.610^^cm-^ (N=9) and n2D = l.TlO^^cm-^ (N=10). Notice 
that the neutral and Hartree bands are shifted with respect to 
each other. The horizontal lines indicates the Fermi energy. 



commodate the extra carriers. However, the electrostatic 
contribution to the gate depends not only on the aver- 
age to density n2D, like in the 2D case, but also on the 
detailed profile qj. As a result, the separation of the 
chemical potential in two contributions, electrostatic and 
quantum mechanical, does not yield a simple procedure 
to 
bo 




FIG. 4: (Color online). Potential ((a),(b)) and density 
((c),(d))profiles, as function of the position across the rib- 
bons, for N = 9 ribbon with n2D = 7.6 10^^ cm~^ (a,c) amd 
n2D ~ 7.6 10^^cm~^ (b,d) The higher density results corre- 
spond to the bands shown in figure 3a. 

In figure 4 we show the self-consistent density profile 
qi and the self-consistent electric potential Ui for a the 
semiconducting ribbon A'^ = 9 with two different n2D ■ We 
see how the density of carriers is larger in the edges than 
in the middle of the ribbon, as expected in a conducting 
system. As a result, the electrostatic potential has a 
inverted U shape. Superposed with these overall trends, 
both qi and Ui feature oscillations, arising from quantum 
mechanical effects. Expectedly, the average Ui is much 
larger in the high density than in the low density. In 
the high density case it is apparent that the potential is 
flat in the inner part of the ribbon, very much like in a 
metal. Very similar trends are obtained for ribbons with 
different widths. 



B. Density vs Vq 

In figure 5a we plot the gate versus the average 2D ex- 
tra density in graphene ribbons with for different widths 
A^, indicated in the figure. This can be a valuable infor- 
mation to estimate the carrier density in gated ribbons 
since Hall measurements in narrow ribbons might be dif- 
ficult to perform. The common feature in all the curves 
is the linear relation between Vq and n2D that reflects 
the dominance of the classical electrostatic contribution 



over quantum effects, exactly like in the 2D case. This is 
also apparent from figure 3, where the shift of the bands 
is much larger than the Fermi energy with respect to the 
bottom of the conduction bands. In addition, we see that 
some of the curves do not intersect at Vg = for n2D = 
in the case of semiconducting ribbons. This is clearly the 
case of = 9 and N = 20. In contrast, the curve TV = 10 
extrapolates to zero. We denote the threshold gate po- 
tential as Vth- and we notice that it corresponds to the 
the change in chemical potential of the ribbon when a 
single electron is added. Since the chemical potential for 
a semiconductor system lies in the middle of the gap and 
the chemical potential when a single electron is added is 
the lowest energy state of the conduction band, we must 
have eVth = where Eg is the gap of the semiconduct- 
ing ribbon. We fit all the curves in figure 5a according 
to: 



n2D{N) = a{N) {VG-Vth{N)) 



(19) 



In figures 5b and 5c we plot a and Vth as a function of 
the ribbon width, N . The wider ribbon considered has 
TV = 60 which corresponds to = 14.5nm. In figure 
5b we see how a rapidly decreases towards the 2D value 
a — —=0.7 10^^cm~^y~^ as W increases. We have ver- 
ified that a scales inversely proportional to the width iV. 
Interestingly, a{N) approaches to the two dimensional 
case very quick, even \i W « d. The evolution of Vth 
as a function of N reflects two facts: on one side, 2 out 
of 3 ribbons are semiconducting. On the other, the gap 
decreases as . The evolution of the gap as a functio: 
of W in the one-orbital tight-binding calculation simila] 
but not identical, to that obtained with density func 
tional calculationsi^. Notice that ribbons oiW = 20nr. 
{N ~ 80) still present a gapi^ and, according to figur 
5b, their capacitance is almost that of two dimension? 
graphene. 



C. Ideal Conductance vs Vb 

We now study the number of occupied bands A/" i: 
a given ribbon as a function of Vq- This is related t 
ideal conductance a function of Vg through the Landaue 
formula for perfectly transmitting ribbons: 



G 



2e' 



(20 



This is equivalent to neglect the effect of disorder in th 
system. Therefore, the value of conductance so obtainC' 
can be considered an upper limit for the real conductanc 
in the system^^. In an ideal armchair ribbon the number 
of chanels Af increases in one by one as Fermi energy with 
respect to the Dirac point is increased. Our approach 
permits to obtain both Af and the quantum shift of the 
chemical potential as a function of the gate voltage and 
to plot A/" as a function of Vq. 

It must be stressed that in the 2D case this disorder 
free model^ accounts for the experimental results'^. This 
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FIG. 5: (Color online), a Gate voltage versus average den- 
sity for armchair ribbons of different widths, b Inverse capac- 
itance, a, as a function of width (measured in terms of A'^). c 
Threshold voltage as a function of the ribbon width A*'. 



might indicate that either disorder is not present in the 
samples or, more likely, because of the suppression of the 
backscattering it does not affect significantly the trans- 
port properties of graphene^. Transport on 2D graphene 
is being extensively studied by a number of groups. The 
effect of imperfections on the transmission of otherwise 
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FIG. 6: (Color online). Conductance as a function of Vg for 
ribbons with = 9 and A^ = 10 (a) and (b), A^ = 40 (c) and 
A' = 60 (d) (see text) 

In figure 6 we plot G(Vg) for various ribbons. In fig. 
6a we plot the conductance G(Vg) for two narrow rib- 
bons {N = 9 and N — 10) in the small gate regime. It 
is apparent that the A^ = 9 ribbon only conducts above 
a threshold gate whereas the A^ = 10 ribbon is conducts 
even for Vg = 0. Therefore, semiconducting ribbons can 
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be electrically tuned from insulating to conducting be- 
haviour with a gate voltage. This is also different from 
the case of 2D graphene which has a rather low resistance 
in the charge neutrality point. In figure 6b we show the 
conductance for the same ribbons at higher Vq- The dif- 
ferent sizes of the plateus for = 9 and iV = 10 reflect 
the different structure of the bands, as seen in figure 3. 
In figures 6c and 6d we show the conductance for wider 
ribbons {N — 40 and iV = 60 respectively). Although 
the shape of the curves is superficially similar for the all 
these ribbons, the Vq necessary to have a fixed number 
of bands at the Fermi energy is a decreasing function of 
W. This is a consequence of quantum confinement: the 
smaller the ribbon the larger the sub-band level spac- 
ing, AEn = E„{k = 0) - £;„_i(fc = 0). A shift of the 
Fermi energy, relative to the self-consistent bands, by an 
amount AEn, so that a new band is available for trans- 
port, implies also an electrostatic overhead, due to the 
change in density, which accounts form most of the eVc, 
as we have discussed above. 

In average, the steps in the stepwise curve A/'(Vg) in- 
crease size as Vg increases. This trend is more apparent 
in wider ribbons. In this sense it can be said that Af is 
sublinear in Vg- In the limit of very wide ribbons, for 
which the density of states is almost 2D, we can derive 
a qualitative relation between number of channels and 
Vq- Since the number bands per unit energy is roughly 
constant, we have that JV scales linearly with the Fermi 
energy with respect to the bottom of the bands. A/" cx Vg. 
On the other side, we know that the density of carriers 
scales linearly with Vg (figs. 2 and 5) and, using the 2D 
density of states we know that density of carriers scales 
with Vg (x n (X Vq. Therefore, we conclude 

Afoiy/Vc (21) 

In figures 6c and 6d we have fitted the numerical data 
to the curve G = aVg, where a and b are fitting pa- 
rameters. We have obtained exponents b — 0.45 and 
b = 0.47 for = 40 and TV = 60 respectively, in agree- 
ment with the qualitative discussion above. In figure 6c 
and 6d we the solid line is the best fit to the equation 
G = a\JVG- Equation ((2T|) is a re-statement of the re- 
lation between Vq and 7^20 of eq. pB]) . In both cases 
the exponent 1/2 comes from the linear relation between 
energy and momentum of the two dimensional carriers 
on one side, and the quadratic relation between n^B and 
Fermi momentum on the other. Therefore, eq. psp and 
pT|) are specific predictions related to the peculiar Dirac- 
like spectrum of electrons in graphene. 

V. DISCUSSION AND CONCLUSIONS 

In this work we study how the electronic structure and 
the density of carriers of graphene based field effect tran- 
sistors evolve as a gate voltage Vg is applied. An impor- 



tant notion is the identification of eVc as the change in 
chemical potential in the graphene layer necessary to ac- 
comodate the density of extra carrers niu- The change in 
chemical potential is the sum of a electrostatic contribu- 
tion, independent of the details of the electronic structure 
of graphene and a quantum contribution that depends on 
the details of the graphene band structure. In the case 
of 2D graphene, we have obtained a particularly simple 
equation (fT8| relating these quantities. We find that the 
quantum contribution is sufficiently small as to go unno- 
ticed in the experiments done so far-i*^"^, but sufficiently 
big as to be measured. Even more important, we claim 
that an independent measurement of uid and Vg would 
provide a direct measurement of the slope of the linear 
energy bands, the so called graphene "speed of light". 

We have considered idealized armchair ribbons which, 
at the charge neutrality point, can be either semicon- 
ducting or metallic. Application of a eVc equal to half 
the band-gap turns semiconducting ribbons into conduc- 
tors. This affords a on/off ratio much larger than that of 
2D graphene. In the case of graphene ribbons the clas- 
sical contribution to capacitance depends on the charge 
density profile qi in the ribbon (fig. 3), which can not 
be measured easily, and we can not provide a simple ex- 
perimental procedure to extract information about the 
electronic structure of the ribbons out of the Vg(?^2d) 
curves. We find that density of extra carriers is higher in 
the edges of the ribbons than in the middle. 

Future work should address the issue of the stability 
of the results obtained in this paper with respect to dis- 
order. In particular, we expect that disorder will turn 
metallic ribbons into semiconducting at Vg = 0. In fig- 
ure 3 it is apparent that the metallic character of the 
ribbon with iV = 10 comes from the states lying in a 
narrow window of width Akx around fc^; = in momen- 
tum states. These states would dissappear if size of the 
ribbon along the x axis is smaller than 

In conclusion, we have studied the electronic structure 
of gated graphene and graphene ribbons in the Hartree 
approximation. This permits to include the Coulomb re- 
pulsion between the extra carriers and their coupling to 
the external gate in a selfconsistent manner. We find 
(eq. (fT8|l a small departure from the classical linear re- 
sult Vg oc n2D which can provide a simple method to 
measure the slope of the graphene bands. In the case of 
semiconducting armchair ribbons we have obtained the 
inhomogeneous distribution of charge carriers qi and po- 
tential Ui along the section of the ribbon. We also find 
that a finite Vg equal to half their band-gap is necessary 
to make them conduct. 
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